Spin-spin interaction and magnetic state of 2-D Wigner crystal 
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It is demonstrated that there is anti-ferromagnetic spin-spin interaction between nearest electrons 
in the 2D Wigner crystal. This is also valid for the Wigner liquid - the state with destroyed long-range 
order but preserved short-range one. We calculate the value of the anti-ferromagnetic interaction 
(both analytically and numerically) and discuss a possible magnetic state of the Wigner crystal. 
This state can be strongly influenced by the spin-Peierls mechanism. 



PACS: 75.10-b, 75.50Ee, 71.10Ay, 73.20Dx, 73.40Qv 



We consider a two dimensional electron gas repulsion) at zero temperature in the presence of a uniform 
neutralizing background. It was shown a long ago by Wigner [Q that at sufhciently low density the electron gas (or 
electron fluid) undergoes a transition into a crystal state. This is because at low density the Coulomb interaction 
dominates the kinetic energy and the correlated state becomes energetically favorable. Analysis of the lattice dynamics 
shows that the stable crystal structure in 2D is the triangular lattice Q . 

The Wigner crystallization has been observed for electrons at the surface of liquid helium |^ . Another 2D system 
for which the electron density can be easily controlled is an inversion layer at a semiconductor surface Q . 

Theoretically the Wigner crystallization has been studied using Monte Carlo simulations, see e.g. Refs. §]. These 
calculations are pretty reliable as far as the critical density is concerned. However there is some controversy about 
possible Ferromagnetic Fermi liquid at a density slightly higher than the crystallization density, see Refs. [0J^. 

Interest in Wigner crystallization has been renewed recently after observation of the insulator-conductor transition 
in dilute 2D electron systems Q . Although this transition probably takes place in the liquid phase it is pretty close 
to the point of crystallization. A very interesting feature of the transition is suppression of the conducting phase by 
in-plane magnetic field pO[ | which influences only spin degrees of freedom. 

In the present work we calculate effective spin-spin interaction in 2D Wigner crystal. This calculation is also valid 
for the Wigner liquid - the state with destroyed long-range order but preserved short-range order. 

To avoid misunderstanding let us note that our calculation does not show any magnetic phase transition in the 
liquid state (i.e. there is no ferromagnetic Fermi liquid between normal Fermi liquid and Wigner crystal). The state 
which we call Wigner liquid is just strongly renormalized normal Fermi liquid. Nevertheless magnetic properties of 
the Wigner liquid (and Wigner crystal) are quite unusual and somewhat similar to that of cuprate superconductors. 
There is competition between superexchange (electron correlation) which gives antiferromagnetic interaction between 
electron spins and the usual exchange Coulomb interaction which gives ferromagnetic contribution. The superexchange 
is proportional to t'^ /U where t is the parameter which describes hopping of an electron to a nearby site and U is 
the Coulomb repulsion for two electrons sitting on the same site. As a result, both the superexchange and exchange 
are proportional to the squared overlap between electron wave functions centered on the different sites of the Wigner 
crystal . Therefore, simple estimates can not answer the question about the sign of spin-spin interaction and we need 
more accurate calculations. To provide better understanding and reliability of the results we have performed these 
calculations twice: analytically and numerically. 

Hamiltonian of the system under consideration is 



1<J 



where pi and are 2D momentum and coordinate respectively. We use effective atomic units which means that all 
distances are measured in units of the effective Bohr radius a*g = 'h^e/m*e'^ , and energies in units of ^s^- Here m* is 
the effective electron mass, and e is the dielectric constant which we assume to be independent of frequency. Number 
density of the electrons n is fixed by condition of clcctroneutrality. An average distance between the electrons 
is defined by nrl = 1/n. It is well established |5| that the crystallization to the triangular lattice ||^ occurs when 
rg Ri 37. In the presence of "disorder" further localization of the electrons stabilizes the Wigner crystal at higher 



densities (rs ~ 10, see g). The distance between the nearest sites in the lattice is equal to a = y 27r/\/3rs sa LQCr^. 
Electrostatic potential acting on the electron near equilibrium position in the lattice is 
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Ui{r)^ const +^^, (2) 
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where r <C a is deviation from the equiHbrium position. To find 7 let us freeze all other electrons in their equilibrium 
positions and calculate Ui{r). Accounting for the six nearest sites gives 7 = 8, and summation over entire lattice 
gives 



aVV-^ 7-9 — r^v? = 5-5171. (3) 



n=l k=0 ^ 

Ground state electron wave function in the potential (0) is 



1 _„2 /o„2 a^/^ 
-1 



i,(r)^^e-^'l^^\ c=^ (4) 



In the above calculation we assume that size of the wave function is much smaller than the lattice spacing, c <C a, 
or (07)^/^ ^ 1. For the crystallization point this parameter equals (07)^/'* — 4. 4, and therefore the approximation is 
well justified in the crystal state. We stress that the parameter appears in the exponent and therefore 4.4 is a very 
large value. Moreover, the approximation is justified in the liquid phase as soon as (07) ^ 1. The matter is that 
the sum (||) is saturated at 2-3 coordination circles and it is independent of the presence or absence of the long-range 
order. For conditions of the experiments (|,|l^ this parameter equals (07)"'^/^ = 3.1. 

To find the magnitude of spin-spin interaction constant (the constant J, which can be substituted to the Heisenberg 
Hamiltonian j) ^i'^j) have to solve a two-particle problem, freezing all the electrons except the nearest two 

ones which are shown by crosses at Figure 1. 
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FIG. 1. Boundary conditions for two-particles problem. 



Hamiltonian of the problem is 

+ f + t/(rx) + C/(r.) + ^, (5) 

where U{y) is potential of all frozen electrons (dots at Fig. 1). The splitting between the ground states for total spin 
S — 1 and S = Q gives us the constant J: 

J^{El=^-El=^) ^{Ea-Es), (6) 

H-^s - Eg-^s. (7) 
H^A - Ea^a- (8) 

Because of the Fermi statistics the two-electron wave function is antisymmetric with respect to permutation. Therefore, 
the symmetric coordinate wave function corresponds to spin S — Q and the antisymmetric one corresponds to 5 = 1. 

Of course the accurate solution of the problem can only be found (and have been found) numerically. However, 
to find the sign and the basic dependence J on a distance parameter a we performed an approximate analytical 
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calculation. To this end, we will follow the procedure, which was suggested long time ago by Gorkov and Pitaevskii 
for the calculation of the term splitting in hydrogen molecule . 

We multiply Eq. (0) by '^a and Eq. (||) by \l/s, take the difference between the results and calculate the integral 
over some region in configuration four dimensional space of the electrons. We choose the integration volume in which 
Xi < X2 (i.e. to the left of the hypcrplanc J^{xi = X2))- Using the Hamiltonian (j^) we obtain 

{Es -Ea) f [ *A*s dri dra = / (^sV*^ - ^aV^-s) dS. (9) 
The kinetic energy term in the right-hand side is reduced to the surface intergral using 

and an integration by parts. 

Now we introduce combinations of the functions ^'1^2 = l/'/SC^s ± '^a)- They correspond to the states of 
"distinguishable" particles, when, e.g. for 5'i(ri,r2), the first electron is principally located near its equilibrium 
position X ~ —a/2 and the second electron near position x = a/2. A simple calculation gives 



*s*A rfri dr2 - i - *2) dv^ dr^ « 1/2. 



Substituting the wave functions 5*1^2 into Eq. and taking into account that under ri ^ y-2 permutation the wave 
functions ^'i ^ 'J'2, we obtain 



J 



*2 

9a; 1 



dx2 dyi dy2. (10) 



The formula ( ]lO| ) shows that the main contribution to the exchange constant is given by the region where the 
electrons are close to each other. Indeed, the x coordinates of both electrons coincide {xi — X2), however, the y 
coordinates may be different. In this case there are strong correlations between the positions of the electrons due to 
Coulomb repulsion. This means that we should go beyond the approximation where the two-particle wave function 
of the electrons is represented as a product of single-particles wave functions. 

It is easy to take into account the effect of the correlations in the quadratic approximation. 

Assuming that the particles are distinguishable and oscillate near their equilibrium positions, we write the Hamil- 
tonian in the following form 

H = -^-^ + Y ((-1 + «/2)' + yl + - al2f + yl) + 

^{|ri-r2| " |ri-a/2| " |r2 + a/2|}- ^^^^ 



Here the frequency u) = ■^/t/o^. 

The Hamiltonian ( p] ) is valid at small displacements Xi and yi from their equilibrium positions. Expanding the 
last term in the curly brackets near (xi,2 — 2^1,2 T a/2) wc finally get the following Hamiltonian in the quadratic 
approximation: 

Ai A2 , UJ^ ( ^2 ,2,-2,2 2_ 



a = - ^ + Y (^^I + yf + ^2 + 2/2 - -ii52 + -2/12/2 ) + 0{xya^). (12) 



Using an obvious change of variables 



2/1 ±2/2 

u, V 



C Xi±X2 , 



V2 ' 
i±. 

we separate Hamiltonian ( p^ into four independent oscillators with frequencies uju,v = -^7(7 ± 1) /a^ and cjc ,, = 
^y {j ^ 2) / . Thus, the ground state wave functions are 
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^2('«, w , ^, ry) = w, ^, -77). 

Substituting Q in Eq. we obtain 

dx2 dy, dy2 = (7 + 2)3/4,/Ia-5/4e-V(7+2)a/2_ 



This formula is presented in the atomic units. In regular units it can be written as 

1/4 



J = 



e 
ea 



a*B 4(7 + 2)^ 



- 1/4 / 








exp ( 


/7 + 2 a \ 


= 3.62— 
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a 



exp 



-1.37, 



(14) 



(15) 



(16) 



The plus sign in the exchange constant shows that the system is anti- ferromagnetic. It is worthwhile to note that the 
exponent in ( p^ ) is different from e~^/^/^, which appears, if states ^'1,2 are represented by a product of independent 
single-particle wave functions (^). 

In order to check the importance of correlations and find the correct exponent for J we have performed numerical 
calculations of the problem over a rectangle area (see Fig. 1). To be absolutely correct we have to impose periodic 
boundary conditions in the rectangle. However, the wave function is very small at the boundary and so the results 
are not sensitive to the boundary condition. It is much more convenient to make the wave function vanish at the 
boundary and this is the condition which we use in the present work. 

The single particle basis set is given by (see Fig. 1) 



(pi{r) = (j)n,n{x',y') 



2 . , 
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(17) 



where A — 3a and B = \/3a. Hence, for the two-electron problem the set is 

\i) = |iii2) = Cilia [0ii(ri)'/'i2(r2)±0ii(i-2)0i2(ri)] 



(18) 



The sign corresponds to S" = (anti- ferromagnetic), and the sign "— " corresponds to 5" = 1 (ferromagnetic). 
The normalization coefficient Ci^i^ = 1/2 if ii = 12 otherwise it equals The matrix element of the Hamiltonian 
(|) is of the for m 



where (j| V('^'2) |j) are matrix elements of the single particle potential and the two-particle interaction correspondingly. 



(1), 



(20) 



To find J, which is exponentially small, we need a very large basis set. The most time consuming part is the 
computation of the two-particle matrix element (j|V(2)|j)j which formally is a 4-dimensional integral. Fortunately 
this integral can be reduced to an integral which is effectively one dimensional. This reduction, which is demonstrated 
in the Appendix, allowed us to perform computations with the size of Hilbert space up to A'^ = 1380. 

Numerical solution of the problem was performed for two different cases. Firstly, we considered all the frozen 
electrons (points at Fig. 1) as point-like charges, which means that the mean-field potential in this case is just the 
sum of the Coulomb potentials: 

Uix,y) = ^uo(|rki - r|), 



uoir) = l/r. 



(21) 
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The sum runs over sites of the triangular lattice. 

Secondly, we considered the density of the frozen electrons to be distributed according to (^) and hence 



Uix,y) = ^ui(|rki - r|), 

kl 

KiM = ^e-'-^/2^^/o(rV2c2), (22) 
c 

where Iq (x) is the modified Bessel function. In both cases all the results are very close and therefore we present plots 
only for the second case. 

The matrices (19) for ferromagnetic and anti- ferromagnetic cases were calculated and diagonalized. The Hilbert 
space was truncated at some high energy state. To be confident that the ground state was found with reasonable 
accuracy we used two basis sets with dimensions N — 975 and N = 1380, where TV is the total number of two-particle 
states. 

The plots of ground state electron density 

p(r) =< 0|,5(r-ri) +,5(r-r2)|0 > (23) 

for a = 45, which corresponds to the limit of our calculations, are given in Figure 2. Similar plots for for a = 15, which 
corresponds to the conditions of the experiments are given in Figure 3. The fact that the maximums coincide 

with the lattice sites tells us about the self-consistency of the method. The shape of the density operator near the 
equilibrium positions also corresponds to the expected Gaussian electron density, obtained from the combinations of 
the wave functions in (p^): 

pix,y)sA = NsA-J, ^'^^"^^"^^ .-^yy' L--.(-+-/2)%e-"-(--'^/2)^ ^2e"'^"'^'/2e-"--'j . (24) 

Here Ns.a = [1 ± e-'^"'^'/2]-i is the normalization coefficient due to the nonorthogonality of the functions ^'i and 
^'2, and u)y = 2w„Wi,/(a)„ -t- Uy), = 2uo^uj^/ {uj(_ -1- uj^). 
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We found that for the whole region of the strength parameter a the anti-ferromagnetic ground state is always below 
the ferromagnetic one. We obtained the following values of the constant J(a) for the parameters of experiments [p|jTc|] 
(e = 8, m* = 0.19rne) 

J(15) = 6.66-10"* = 0.6K, 
J(45) = 1.48 • 10"*^ 1.4 • IQ-^K. 

(25) 

Experiments HJl^ correspond to a ~ 15. The behavior of J, as expected, has an exponential dependence ~ e~^^. 
The plot of the dependence ln(J) vs. ^/a summarizes our results in Figure 4. Diamonds and crosses show the 
magnitude for different basis sets and nicely depict the truncation effects for large a and for small number of basis 
states. The dot-dashed line represents the theoretical curve (|l^) and the solid line is the best fit. 




FIG. 4. Dependence ln(J) on ^/a. Diamonds and crosses represent data for large (A'' = 1380) and shortened (A'^ = 975) 
bases correspondingly. The dot-dashed line shows the theoretical curve of Eq. (|l5|). The solid line is the best fit by curve 
y = C - Sy/k- l3ln{a). 
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We fitted our data by two different functions. For the first one we fixed the power of a in the preexponential factor 
(the analytical formula ([l5| ) gives a~^/^): 

ln( J) = C - S^/^ - 5/2 In(V^) 

and got C = 2.25, S — 1.595. Let us note that the constant 6 is slightly larger than the predicted analytical value 
San = V7 + 2/2 = 1.37. This fact shows that the higher order terms in the expansion of the Hamiltonian (jl^) are 
important for getting the correct magnitude of the exponent. 
In the second case we looked for the best parameters for 

\n(J) ^ C - - (3\n{y^). 

We obtained C = 2.02, 5 = 1.74, /3 = 1.88. 

The analytical formula for J in Eq. ( p^ ) and more accurate fits of the numerical calculation data allows one to 
estimate the value of J for the Wigner crystal and Wigner liquid states in the large region of densities. 

Due to geometric frustration coUinear long-range antiferromagnetic order is, strictly speaking, not possible on a 
triangular lattice. The possible solution in this case is a system spin wave function which in the zero approximation 
consists of spin zero pairs. The antiferromagnetic interaction increases when the distance between the electrons 
decreases. Therefore, there should be a tendency for nearby electrons coupled to spin zero pairs to move slightly 
toward each other. This phenomenon is usually called the "Spin-Peierls" mechanism. 

It is interesting that the value of the spin-spin interaction J — J(15) is comparable with the energy in the 
experiments |^,|l3, where H is the critical magnetic field destroying conductivity. We can speculate that this field 
effectively transforms the system to a ferromagnetic state. In the ferromagnetic state the conductivity should be 
smaller than in the antiferromagnetic state. Indeed, hopping of the electrons with opposite spins from one site to 
another is allowed by the Pauli principle. The magnetic field rearranges spins in the same direction. In this case such 
hopping is suppressed by Pauli blocking. This possibly destroys the conductivity. 
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APPENDIX A: MATRIX ELEMENTS OF THE INTERACTION 



The eigenfunctions and eigenvalues for the single-electron problem in the two-dimensional rectangle with sides 
A = Sa and B = ^/Sa are given by Eq. (p^. The matrix element of the Coulomb interaction between two unfrozen 
electrons is given by the integral over the rectangle's area 

ViM,^ f <PMY<PM) \, \, M r'i)^M)d^r[d\2 (Al) 



It is convenient to change variables (see Fig. 1): 

II I I II I II I / A o\ 

u = x^-X2, v^y^-y2, X2 = X2, 2/2=2/2 (A2) 
Introducing notations u = ^ , v — ^ , and 

-f^nin2n3n4(ii,52) = 4 sin(£ini) sin(xin3) sin(i2'^2) sin(i2'^4), 

nTT — U 

^^711 712^3714 (^) ~ / ^711712 713 714 (*^2 ~\~ U., X2^dX2 



Jn.n.nsn, = ^ (l + (-1)"^+"^+"^+"^) , (A3) 

the matrix element ( |Al[ ) can be rewritten as: 

A _ _ dudv 

ABtt^ Jo Jo V u + w 
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In this transformation we use the following relation 



(X, X + U)dx= (_l)"i+"2+n3+ru 



nin2n3"4 



(x + u, x) dx 



It is convenient to calculate the double integral (A4) using polar coordinates u — ^ — -^cosi, v = ^ — rsint. 
Taking into account that tan^^{B /A) = 7r/6 we find that 



V, 



V2 = 



7r/6 pTT\/Z/ cos(t) 



(it 



r/3 



(it 



tt/ cos(i) 



r cos{t) / \/i)Wrmra2ni3rai{r sin(t)) dr, 
(r sin(t)/\/3) Wmim2„3„4 (r cos(t)) dr. 



(A5) 



In order to write down an analytical expression for the function Wnin2n'ini{u) let us introduce the following notations 



Using 



one can find that 



^^nin2n3n4 ('^) 



n ~ \n'i — ni\, n = 0;l,2,. 

m = |n4 — ri2|, m = 0,1,2,. 

^ = ris + rti, ; = 2, 3, . . . 

fc = n4 + rt2, fc = 2, 3, . . . 



f2(n,l,k]u) if n = 7Ti ^ 

/2(Z, n, m;u) if I = k 

— f2(n,m,l;u) \i n — k 

—f2{m,n,k;u) if m — I 

f3(l,k;u) if n — m — 

fi{n,l;u) if n = m =/= and I = k 

fbO'j if n = TO = and / — k. 



(A6) 



= < 



(A7) 



where 

fi{n,m,l,k;u) 
f2{n,m,l;u) 
f'i{n,m;u) 
fi{n,m;u) 



sin(nw) (/c^ — to^) m sin(MTO) (/^ — n?) I sin(Zu) (A;^ — to^) fc sin(ufc) (^^ — n^) 



{in? — fc2) (n^ — m?) {m? — P) {m? — n^) (/2 _ ^2") _ ^2^ (— + k"^) [k"^ — P) 

— ■m?P \ msin(MTO) (n^ — Z sin(Zu) (n^ — m^) 



TT — u sin(nM) 
— - — cos(nM) + 



2n 



1 + 2- 



(n? — m?) [m? — P) J {irP — P) [m? — rP) [P — rP) {P — rrP) 



= TT — M — 



sin(uTO)n^ sin(riw)TO^ 



[nP — rP) TO (tt.^ — TO^) n 



TT — u sinfnu) 3rP + rrP sin(uTO) 3 irP + rP 
[cos(nu) + cos(uto)J ' ' 



In — m-^ 



2m — 



fz{n]u) = TT — M H ^ — cos(n7i) + — sin(nu) 



(A8) 



This completes the description of the calculation procedure for the two-electron Coulomb matrix element. The 
advantage is that each of the two integrals in (A5) require numerical work equivalent only to the computation of a 
ID integral. 

Calculation of the the single-particle matrix clement of the external potential U{v) is much simpler. It is convenient 
to use x and y instead of x' and y' (see Fig. 1). Then 



11 13 



{nimi, 7137713} J^j^ 



4 

— — (— 1) 2 JninsJmims 



A/2 /•B/2 



U{x,y) X 



cos(— nx) - (--1)"" cos( — |cos(— my) - (-1)™=^ cos{ — ky)jdxdy, 



(A9) 
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where 



n 
m 
I 

{ k 



\n3 - nil 
|m3 - mil 
ns + ni 
ms + mi 



(AlO) 
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